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At the filling factor v=2, the bilayer quantum Hall system has three phases, the spin triplet phase, 
the spin singlet phase and the canted antiferromagnet (CAF) phase, depending on the relative 
strength between the Zeeman energy and interlayer tunneling energy. We present a systematic 
method to derive the effective Hamiltonian for the Goldstone modes in these three phases. We then 
investigate the dispersion relations and the coherence lengths of the Goldstone modes. To explore 
a possible emergence of the interlayer phase coherence, we analyze the dispersion relations in the 
tunneling energy zero limit. We find one gapless mode with the linear dispersion relation in the 
CAF phase. 
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I. INTRODUCTION 

In the bilayer quantum Hall (QH) system, a rich 
physics emerges by the interplay between the spin and 
the layer (pseudospin) degrees of freedom[l|, For in- 
stance, at the filling factor v = 1, there arises uniquely 
the spin-ferromagnet and pseudospin-ferromagnet phase, 
showing various intralayer and interlayer coherent phe- 
nomena. On the other hand, the phases arising at v = 2 
are quite nontrivial. According to the one-body pic- 
ture we expect to have two phases depending on the 
relative strength between the Zeeman gap Az and the 
tunneling gap Asas- One is the spin-ferromagnet and 
pseudospin-singlet phase (abridged as the spin phase) 
for Az > Asas; the other is the spin-singlet and pseu- 
dospin ferromagnet phase (abridged as the ppin phase) 
for Asas > Az. Instead, an intermediate phase, a canted 
antiferromagnetic phase (abridged as the CAF phase) 
emerges. This is a novel phase where the spin direc- 
tion is canted and make antiferromagnetic correlations 
between the two layers [3, |4j. Das Sarma et al. obtained 
the phase diagram in the Asas — d plane based on time- 
dependent Hartree-Fock analysis, where d is the layer 
separation^ Later on, an effective spin theory, a 
Hartree-Fock-Bogoliubov approximation and an exact di- 
agonalization study were employed to improve the phase 
diagram [5|-|8j]. Effects of the density imbalance on the 
CAF were also discussed [gL [Toj[ . 

The first experimental indication of the CAF phase was 
given by inelastic light scattering spectroscopy [ll|. They 
also have observed softening signals indicating second- 
order phase transitions [12]. Subsequently, an unambigu- 
ous evidence of the CAF phase was obtained through 
capacitance spectroscopy as well as magnetotransport 
measurements p^j - [T§ ] . 

The ground state structure of the v = 2 bi- 
layer QH system has been investigated based on the 
SU(4) formalism 20- 25]. The expectation values of the 
SU(4) isospin operators are the order parameters, in 
terms of which an anisotropic SU(4) nonlinear sigma 



model has been derived to describe low-energy coherent 
phenomena[20]. However, the effective Hamiltonian for 
the Goldstone modes has not been derived. Though there 
are some results with the use of Grassmannian fields in 
the spin and ppin phases, no attempts have been made 
in the CAF phase. On the other hand, experimentally, a 
role of a Goldstone mode has been suggested by nuclear 
magnetic resonance [19( in the CAF phase. 

In this paper we develop a generic formalism to deter- 
mine the symmetry breaking pattern and to derive the ef- 
fective Hamiltonian for the Goldstone modes in the three 
phases of the v = 2 bilayer QH system. The symmetry 
breaking pattern reads 



SU(4) U(l) ® SU(2) ® SU(2), 



(1.1) 



and there appear eight Goldstone modes in each phase. 
The corresponding Goldstone modes in the two phases 
match smoothly at the phase boundary. All the modes 
are actually gapped except along the phase boundaries 
due to explicit symmetry breaking terms. It is im- 
portant if gapless modes emerge in the limit Az — > 
or Asas — > 0, where the spin coherence or the inter- 
layer coherence is enhanced. Gapless modes are genuine 
Goldstone modes associated with spontaneous symmetry 
breaking. Naturally we have gapless modes in the spin 
phase as Az — > and in the ppin phase as Asas 0. 
It is intriguing that we find one gapless mode with the 
linear dispersion relation in the CAF phase as Asas — > 0. 

This paper is organized as follows. In Sec. UH we re- 
view the Coulomb interaction of the bilayer QH system 
projected to the lowest Landau level (LLL) and the SU(4) 
effective Hamiltonian after making the derivative expan- 
sion. We also review the ground state structure in the 
three phases. In Sec. IIII1 which is the main part of this 
paper, we develop a unified formalism to derive the ef- 
fective Hamiltonian for the Goldstone modes. Then we 
discuss the SU(4) symmetry breaking pattern and the 
Goldstone mode spectrum, such as the dispersion rela- 
tions and the coherence length in each phase. In partic- 
ular, for the investigation of the CAF phase, we find it 
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useful to introduce two convenient coordinates of SU(4) 
group space, the s-coordinate and the p-coordinate. We 
study the dispersions and the coherence length in the 
limit Asas — > 0, to explore a possible emergence of the 
interlayer phase coherence in the CAF phase. Remark- 
ably, we find one coherent mode whose coherence length 
diverges. Section ITVl is devoted to discussion. 



II. THE SU(4) EFFECTIVE HAMILTONIAN 
AND THE GROUND STATE STRUCTURE 

Electrons in a plane perform cyclotron motion under 
perpendicular magnetic field B± and create Landau lev- 
els. The number of flux quanta passing through the sys- 
tem is N<s> = B±S/$n, where S is the area of the sys- 
tem and $d = 2nh/e is the flux quantum. There are 
iV<i> Landau sites per one Landau level, each of which is 
associated with one flux quantum and occupies an area 
S/N& = 2n£ 2 3 , with the magnetic length £b = yh/eB±. 

In the bilayer system an electron has two types of in- 
dices, the spin index (f,i) and the layer index (f, b). 
They can be incorporated in 4 types of isospin index a — 
ft,£J.,bt,bJ.. One Landau site may contain four electrons. 
The filling factor is v = N/N$ with N the total number 
of electrons. 

We explore the physics of electrons confined to the low- 
est Landau level, where the electron position is specified 
solely by the guiding center X = (X,Y), whose X and 
Y components are noncommutative, 



-U%. 



[X,Y] 

By introducing the ladder operators, 



(2.1) 



b =vk< x -' n b ' = ^F a (x + in (2 - 2) 

obeying [6, b'] = 1, we construct the Fock states, 

L (&t)»|0), n = 0,1,2,..-, 6|0)=0. (2.3) 



These states are the Landau sites in the symmetric gauge. 

We expand the electron field operator by a complete 
set of one-body wave functions ip n (x) = (x\n) in the 
LLL, 



ll> a {x) 



E 

71=1 



c a (n)ip n (x), 



(2.4) 



where c a (n) is the annihilation operator at the Landau 
site |n) with a = ft,£J,,bt,bi. The operators c Q (m), ct(n) 
satisfy the standard anticommutation relations, 

{c a (m), c^(n)} = 5 mn 5 a p, 

{c a (m),cp(n)} = {4(m),4(n)} = 0. (2.5) 



The electron field ip a (x) has four components, and the 
bilayer system possesses the underlying algebra SU(4) 
with having the subalgebra SU sp i n (2) x SU pp i n (2). We 
denote the three generators of the SU sp i n (2) by r„ pm , 
and those of SU pp i n (2) by rP pm . There are remaining 
nine generators T^ pin T^ pm . Their explicit form is given in 
Apendix A. 

All the physical operators required for the description 
of the system are constructed as the bilinear combina- 
tions of tp( x ) and 4>^(x). They are 16 density operators 



p(x) = 


i/;\x)ip(x), 


S a (x) = 


^(x)r^(x), 


Pa{x) = 


^(x)rryj(x), 


Rab(x) = 





(2.6) 

where S a describes the total spin, 2P Z measures the 
electron-density difference between the two layers. The 
operator i? Q f, transforms as a spin under SU sp i n (2) and 
as a pseudospin under SU pp i n (2). 

The kinetic Hamiltonian is quenched, since the kinetic 
energy is common to all states in the LLL. The Coulomb 
interaction is decomposed into the SU(4)-invariant and 
SU(4)-noninvariant terms 

H c=lJ d 2 xd 2 yV+(x - y)p(x)p(y), (2.7) 

H c= 2 J d 2 xd 2 yV-(x - v)P z {x)P z (y), (2.8) 



where 



V ± (x) 




(2.9) 



with the layer separation d. The tunneling and bias terms 
are summarized into the pseudo-Zeeman term. Combin- 
ing the Zeeman and pseudo-Zeeman terms we have 

ffzpz = - J d 2 x(A z S z + A SAS P X + A bias P 2 ), (2.10) 

with the Zeeman gap Az, the tunneling gap Asas, and 
the bias voltage A bias = eVb ias - 
The total Hamiltonian is 

H = H+ + H c + H ZpZ . (2.11) 

We investigate the regime where the SU(4) invariant 
Coulomb term Hq dominates all other interactions. Note 
that the SU(4)-noninvariant terms vanish in the limit 
d, A z , A SAS , A bias -> 0. 

We project the density operators (|2.6p to the LLL by 
substituting the field operator (|2.4p into them. A typical 
density operator reads 



Rab{p) 



1 *-a< 



(2.12) 
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in the momentum space, with 

^(p) = ^E("i e " pX i m ) ct (") T « pinr r in c(m), 

rrm 

(2.13) 

where c(m) is the 4-component vector made of the oper- 
ators c a {m). 

What are observed experimentally are the classical 
densities, which are expectation values such as p cl (p) = 
(S\p(p)\S), where \S) represents a generic state in the 
LLL. The Coulomb Hamiltonian governing the classical 
densities are given bv(23j 



# eff = 7T J d 2 pV+(p)p C \-p) P C \p) 

+ 47T / d 2 pV D (p)Pf(-p)Pf(p) 



d' P vz(p)[s: i (-p)s c a i (p) + p c a \-p)P c a \ P ) 



+ Rf b (-p)Rf b (p)]-n J d 2 pV x (p)[S:\-p)S:\p) 

+ Pf(-p)Pf(p) + Rt{-p)RtiP)] 

-\J ' d 2 pV x (p)p*(-p)p c \p), (2.14) 



Hamiltonian density 

n cS = jf (^T(d k s a ) 2 + (d k Va) 2 + (d k n a , 
+ 2j- (J2( d kSa) 2 + (d k V z ) 2 + (d k TZ az ) 2 ) 

+ P4>[ec ap (V z ) 2 - 2e x (]T(^) 2 + (K az ) 2 ) 

-(4 {Sa) 2 + {r a ) 2 + {n ab ) 2 ) 

- (A Z S Z + A SAS V X + A Kms V z ) - (4 - e x )}, 

(2.17) 

where p$ = p^jv is the density of states, and 
1 



J.s 



16V27T 



jd _ J 



2 d 

IT in 



1 + ^ 2 /2i% eTic ( d /^ B \ 



tx 



1 Fk 



2 V 2 



1 ± e «P/a&erfc (d/v^a) 



— -T7T~ ^Gi £ cap — 4e n — 2e x , 



(2.18) 



with 



where Vd and Fx are the direct and exchange Coulomb 
potentials, respectively, 



Vd(p) 
Vx(p) 



-e 2 BP 2 /2 



4ire\p\ 
V2^e 2 £ B 

47T£ 



I (£ 2 B p 2 /4)e- e > 2 /\ (2.15) 



and 



= w ( 



V?(p) = 
± 



27re 2 £ B 
8^ 

e 2 £ 2 ,oo 



J (4p 2 /4)e-^ p2/4 



47T6 



i; 1 dke-^ k2 ~ kd J (i 2 B \p\k), (2.16) 



with Iq(x) the modified Bessel function, and Jo{x) the 
Bessel function of the first kind. 

Since the exchange interaction V ± (p) is short ranged, 
it is a good approximation to make the derivative expan- 
sion, or equivalently, the momentum expansion. We may 
set p c \p) = pa, S*(p) = p$5 a (p), P a cl (p) = pMp), 
and R a \(p) = p^'R-abip) for the study of Goldstone 
modes. Taking the nontrivial lowest order terms in 
the derivative expansion, we obtain the SU(4) effective 



(2.19) 



This Hamiltonian is valid at v = 1 and 2. 

It is to be remarked that all potential terms vanish 
in the SU(4) invariant limit, where perturbative excita- 
tions are gapless. They are the Goldstone modes associ- 
ated with spontaneous breaking of the SU(4) symmetry. 
There are eight Goldstone modes, as we shall show in 
Section [TTT1 They get gapped in the actual system, since 
the SU(4) symmetry is explicitly broken. Nevertheless 
we call them the Goldstone modes. 

The ground state is obtained by minimizing the ef- 
fective Hamiltonian (|2.17[) for homogeneous configura- 
tions of the classical densities. The order parameters are 
the classical densities for the ground state. It has been 
shown[20] at v = 2 that they are given in terms of two 
parameters a and /3 as 



As 



V 
111 



A S AS 2 

a 



Ac 

Asas 



o _ A S as „ 2 



A 



VP, 



as/1 - a 2 /3, 



Kv = -^Wi-^Vw 2 , 



(2.20) 
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with all others being zero. The parameters a and /3, 
satisfying |a| < 1 and \/3\ < 1, are determined by the 
variational equations as 



A 2 = A 2 AS 4 £x (A 2 -/3 2 A 2 AS ) 



1 - /3 2 Ao^T^W 
Abias _ 4 (ejg + 2a 2 (e D - e x )) 1 



A 



VW 1 



(2.21) 
= , (2.22) 



/3A SAS 
where 



A = ^A 2 AS a 2 + Al(l - a2)(l - /3 2 ). (2.23) 

As a physical variable it is more convenient to use the 
imbalance parameter defined by 



_ ^,0 A S AS 2o 



(2.24) 



instead of the bias voltage Abi as - This is possible in the 
ppin and CAF phases. The bilayer system is balanced 
at (To = 0, while all electrons are in the front layer at 
(T = 1, and in the back layer at a = — 1. 

There are three phases in the bilayer QH system at 
v = 2. We discuss them in terms of a and (3. 

First, when a = 0, it follows that 5° = 1, V° = K° ab = 
0, since A = A Z \/l — (3 2 . Note that (3 disappears from 
all formulas in (|2.20p . This is the spin phase, which is 
characterized by the fact that the isospin is fully polar- 
ized into the spin direction with 



1, 



(2.25) 



and all others being zero. The spins in both layers point 
to the positive z axis due to the Zeeman effect. 

Second, when a = 1, it follows that S° z = and ("P°) 2 + 
(V®) 2 = 1. This is the ppin phase, which is characterized 
by the fact that the isospin is fully polarized into the 
pseudospin direction with 



V° x = V^P, Vl = (3 



Co, 



(2.26) 



and all the others being zero. 

For intermediate values of a (0 < a < 1), not only 
the spin and pseudospin but also some components of the 
residual spin are nonvanishing, where we may control the 
density imbalance by applying a bias voltage as in the 
ppin phase. It follows from (|2.20l) that, as the system 
goes away from the spin phase (a = 0), the spins begin to 
cant coherently and make antiferromagnetic correlations 
between the two layers. Hence it is called the canted 
antiferromagnetic phase. 

The interlayer phase coherence is an intriguing phe- 
nomenon in the bilayer QH system [lj]. Since it is en- 
hanced in the limit AgAs — > 0, it is worthwhile to inves- 
tigate the effective Hamiltonian in this limit. We need to 
know how the parameters a and (3 are expressed in terms 
of the physical variables. Form (12.21)) it is trivial to see 



that (1 - /3 2 )/A§ AS = 0(1). Up to the order 0(A| AS ), 
([2~2"2"|) is reduced to 



l-(3< 



The solutions are 



A!(i-« 2 ) + %^ 2 



= 0. 

(2.27) 



A ^ ] ' 0(A SAS ), (2.28) 



with 



A ^ A SAS + 0(A3 Ag ), (2.29) 
for ([2~2U1) . By using (j2~2"4]l we have 



V: 



(Tq = ±a 




0(A 4 S 



SASJ 



(2.30) 



The parameters a and /3 are simple functions of the phys- 
ical variables Ag A s/Az and <To in the limit Ag A s — > 0. 

In particular, one of the layers becomes empty in the 
ppin phase and also near the ppin-phase boundary in the 
CAF phase, since we have cto — > ±1 as a — » 1. On the 
other hand, the bilayer system becomes balanced, since 
we have uq — > as a — ¥ in the spin phase and also near 
the spin-phase boundary in the CAF phase. We might 
expect novel phenomena associated with the interlayer 
phase coherence in the CAF phase. 



III. EFFECTIVE HAMILTONIAN FOR 
GOLDSTONE MODES 

Having reviewed the three phases in the bilayer system 
at v = 2, we proceed to discuss the symmetry breaking 
pattern and construct the effective Hamiltonian for the 
Goldstonc modes in each phase. There is a systematic 
method for this pur pos e, which was developed in particle 
and nuclear physics [261 \27\. 

We analyze excitations around the classical ground 
state (|2.20p . It is convenient to introduce the SU(4) 
isospin notation such that 



AO) 

-a0 



oO T (0) _ .pO T (0) _ TpO 



(3.1) 



(0) 



We set all of them into one 15-dimensional vector X^J 

with the index \xv: Note that there is no component Iqq . 
Most general excitations are described by the operator 



C\p | I ^ TT^sTjS 
76 



Z°>„>, (3-2) 



where T 7 a are the matrices of the broken SU(4) genera- 
tors in the adjoint representation of SU(4), each of which 
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is a 15 x 15 matrix. The greek indices run over 0, x, y, z. 
The phase field tt 1 s{x) are the Goldstone modes associ- 
ated with the broken generators, and the coefficient X{x) 
is the amplitude function corresponding to the "sigma" 
field in the linear sigma model. 

It has been argued (2p| that there are nine indepen- 
dent real physical fields. They are the amplitude fluctu- 
ation field X(x) satisfying X 2 (x) < 1, and eight Goldstone 
modes ■k 1 s(x). Hence, only eight generating matrices T-ys 
are involved in the formula (|3.2p . We shall explicitly de- 
termine them in each phase in the following subsections. 
Since we are only interested in an effective low energy 
theory of the Goldstone bosons, we set X(x) = 1. Then 
we may identify 

S a =X a Q, ~P a = loa, TZ-ab — X a b, (3-3) 

and express various physical variables in terms of the 
Goldstone modes 7r 7 a(x). 

We expand the formula (|3.2j) in tt^s, 

V (a;) = + Z# (x) + 1$ (!) + ■••, (3.4) 

where xffi (x) is the nth order term in the Goldstone 
mode T:-ys- Up to the second order, they are 

-^jj,J i x ) = ^| fnv^S,^ 1 !/" f/j,"v" ,-y'6' ,tai'Ky6' ! {'Y'5'3?jl u , (3-6) 

where f^jS^'v' are the structure constant of SU(4), 



(3.7) 



about which we explain in Appendix A (|A7[) . 

Each phase is characterized by the order parameter 
, which are nothing but ([2.20j) . The key observation is 

that the first order term (x) contains all informations 
about the symmetry breaking pattern and the associated 
Goldstone modes, yielding their kinematic terms. On the 
other hand, the second order term X^J (x) provides them 
with gaps. 



A. Spin Phase 

First we analyze the spin phase. Setting a = in the 
order parameters (|2.20l) . we obtain 



(3.8) 



With the use of this, it is straighforward to calculate the 
first order term XpJ (x) in (13.51) , 



X^ - -it l( 



-'.Ijl 



1 XfJ,- 



(3.9) 



There are eight fields ir yil and ir xfl with \i = 0, x, y and 
z, which are the Goldstone modes. Since they emerge 



in eight directions, x/i and yu, the broken generators 
are T xjl and T y/Jt . Consequently, the symmetry breaking 
pattern reads 

SU(4) -> U(l) Za0 <» SU(2) Ioa ® SU(2)x 2a , (3.10) 

implying that the unbroken generators are T z q, To a and 
T 

± za ■ 

We require p. 21) to satisfy the SU(4) algebraic relation 

[X afJ ,(x,t),Xb IJ ,(y,t)] = ie abc p^ 1 X c0 (x,t)S{x - y), (3.11) 

so that the field X^ v describes the SU(4) isospin. From 
p. lip , we obtain the equal-time commutation relations 
for the Goldstone modes, 



Tr Xf t(x, t), n vfl (y, t)] = iS(x - y), 



(3.12) 



with 7r 7 5 = tt^s- Equivalently, we may construct a 
Lagrangian formalism so that (|3. 12|) is the canonical com- 
mutation relation. 

It follows from (|3.3p and (|3.9p that the eight Goldstone 
modes are explicitly given by 

$x TTyO; ^xQ j ^xa TTya; T^ya — ^xa- 

(3.13) 

Substituting them into (I2.17p . we obtain the effective 
Hamiltonian of the Goldstone modes in terms of the 
canonical sets of tt x ^ and n yil as 



2J. 



P0 

Po 



fl— 0,2 



^2 [(dkTTxa) 2 + {dk^yaf 



a—x,y 
-(2) 



(3.14) 



(2.) ______ 

where IqJ are given by (|3.6p , and read 



T K ' — 



fi=0,x,y,z 



r (2) 'Kxx'^yy ^yy^xx ^xy^y. 



x n yx n xy 



'-0z 



(3.15) 



The annihilation operators are defined by 



m = 
m = 



TT x0 + ITTyQ 



V2 ' V2 ' 

and satisfy the commutation relations, 

Vi{x,t),r]Uy,t) =5ijS(x-y), 



i yx 



V2 
tt X7 + in,. 



(3.16) 



(3.17) 
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with i,j = 1,2,3,4. 

The effective Hamiltonian (|3.14p reads in terms of the 
creation and annihilation variables (|3.16[) as 

H s PiD = ^y d kV id kVa + ^ V d k iid kVa 



Po 



Po 



a=2,3 



A Z ^ Vila 



a=l,4 



[A 



A 



[^2% - vim] 



a=2,3 

A 



—[4 m - vim,]- 

(3.18) 



The variables 772, 773 and 774 are mixing. 

In the momentum space the annihilation and creation 
operators are 77^ and ri\ k together with the commuta- 
tion relations, 



5ij8(k — fc') 



(3.19) 



For the sake of the simplicity we consider the balanced 
configuration with Abias = in the rest of this subsec- 
tion. Then the Hamiltonian density is given by 



d 2 k H sp " 



where 



-T/spin 

H 2 



/T/spin 
ft* 



^k 2 

Po 

A T d 
-^k 2 

Po 

A T d 
-^k 2 

Po 



Ac 



A 



Ac 



x 



4J, 



A 



SAS 



(3.20) 

(3.21) 
(3.22) 
2 + A z 
(3.23) 



We first analyze the dispersion relation and the coher- 
ence length of 771 ,fe. From (|3.2ip . we have 



E m (k) 



4J. 
Po 



s k 2 



(3.24) 
(3.25) 



The coherent length diverges in the limit Az — > 0. This 
mode is a pure spin wave since it describes the fluctuation 
of S x and S y as in p.!3|) . Indeed, the energy Q3.24p as 
well as the coherent length p. 251) depend only on the 
Zeeman gap Az and the intralayer stiffness J s . 
We next analyze those of 772, fc, 



E m (k) 



*V</2 



A J d 
^k 2 
Po 

2Jjm/t 



Ac 



X ! 



Ac 



(3.26) 
(3.27) 



A' 



They depend not only Az but also on the exchange 
Coulomb energy c x and the interlayer stiffness originat- 
ing in the interlayer Coulomb interaction. 

We finally analyze those of 773^ and 774^, which are 
coupled. Hamiltonian (|3 . 231) can be written in the matrix 
form, 



/T/spin 
H 3 



m,k 



A k -iAsAs 
«A SA s B k 



V3,k 



where 

4 T d 
A k = ^k 2 
Po 



Ac 



x> 



B k = ^k 2 + A z 
Po 



Hamiltonian (|3.28p can be diagonalized as 

t 

ain _ / 77 3 \ 

L 3 



-T/spin 

tto — 



m,k 
m,k 





£ f ' 4 



(3.28) 
(3.29) 

(3.30) 



where 
E^ = 

E^ = 



(A k - B k ) 2 + AA 2 



SAS 



A k 



B k — \l {A k — B k ) 2 + 4A| AS 



(3.31) 



and the annihilation operator 77^ (7 = 3, 4) given by the 
form 



1 I y/C'i + AA 2 



SAS 



Cfc ) 773^fc - 2A S AS?74,fc 



h (C 2 + AA 2 AS + C fev ^+4A|^) 
-* (VCfc + 4A^ AS - Cfc) 773^ + 2A SA s77 4!fe 



2 C2 + 4A§ AS -<W<7ji + 4Ag AS 



(3.32) 



with C k = A k — B k . The annihilation operators (|3.32p 
satisfy the commutation relations 



Vi,k,fjlk' 



S^Sik-k'), 



(3.33) 



with 7, j = 3,4. We obtain the dispersions for the modes 
fj i<k (i = 3, 4) from (|3~29j) and (j33T]) 

By taking the limit k — ► in (|3.31l) . we have two gaps 



^k^O 
r"J4 _ 



A z + 2c x + [A(e- X ) 



- \2 
2 



Asas] 



(3.34) 



= + 2e x - [4(e^ + Ag AS J 2 . 

The gapless condition (E^o = 0) implies 

A z (Az+4e x )-Ai AS = 0, (3.35) 

which holds only along the boundary of the spin and 
CAF phases: See (4.17) in Ref.[20|. In the interior of 
the spin phase we have Az(Az + 4e A ) — Ag AS > 0, as 
implies that there arise no gapless modes from 773 and 7)4. 
These excitation modes are residual spin waves coupled 
with the layer degree of freedom. 



7 



B. Ppin Phase 

We next analyze the ppin phase. Setting a — 1 in the 
order parameters (|2.20|) . we obtain 



(3.36) 



In order to determine the symmetry breaking pattern, 
we rotate this vector around the Oy axis so that only one 
component becomes nonzero. We can show that 



by choosing 

V p (0f,) = eMiOpT Qy ), 



(3.37) 
(3.38) 



with cos dp = yl — (3 2 and s'mdp — —(3. 

In the rotated basis the order parameter has a single 
nonzero component just as (|3.8p in the case of the spin 
phase. Therefore the further analysis goes in parallel 
with that given in the previous subsection. Namely, there 
are eight Goldstone fields, 



TP(1) - 



flZ 5 



JP(1) = n P 



and the symmetry breaking pattern reads 
SU(4) ->■ SU(2) X P n ® U(1) x p ® SU(2) x p 



(3.39) 
(3.40) 



precisely as in the spin phase. 

Let us relate the variables in the rotated system to the 
original variables in the formula (|3.2I) . The SU(4) isospin 
operator after the rotation is given by 



W = [^(9CV(4 



(3.41) 



with the use of (|3.38j) . We substitute (|3.2j) into this for- 
mula to find 



7<5 



■p(0) 
li'v 1 ' 



with (|3.37|) . where 7rJ? a is defined by 



rP(0) 
L 7 <5 



[Wis)] 



7'*' tO 

7 5 -^r'i'' 



(3.42) 



(3.43) 
(3.44) 



while l^g has been used by (|3.37|) . Here, we have used 
the formula of the SU(N) group, 

TbK =J2 Tb [exp*0.Ad(T o )]£ $ c 

b b 

= exp [iO a T a ] $ 6 T b exp [-i6 a T a ] , (3.45) 

where $b is an arbitrary adjoint vector with a, b, c — 
1, . . . , dim SU(N), and exp [i9 a T a ] is the element of 
SU(N). Here we have N= 4 and corresponds to tt^. 



The SU(4) isospin density fields l p v satisfy the SU(4) 
algebraic relations 



x la{ x ,t)^ x lb{y^) =ieabcP<t, 1 Zo c (x,t)5(x-y), (3.46) 



from which we obtain the canonical commutation rela- 
tions for the Goldstone modes, 



x-y) 



(3.47) 



with n» v = py 2 ^. 

We go on to derive the effective Hamiltonian governing 
these Goldstone modes. The first step is to convert the 
relation (|3.4ip to express the original fields in terms of 
those in the rotated system. Explicitly we have 

I a0 =I a p , l w =^. (3.48) 
The second step is to expand f|3 .42[) in terms of ir^ s , 



V = 
m 


-^ 2 +0(tt 2 ), 2*, = **,+ 


0(tt 2 ), 




7T P 7T P + 

zz yy ' 


7T P 7T P - 

yy zz 


7T P 7T P - 

zy n yz 






2 






JP _ 

x y0 ~ 


7T P 7T P + 

zy xz 1 


7T P 7T P - 

xz zy 


7T p 7T p - 

zz xy 


7T P 7T P 

2:2: 




2 






7 p _ 


7T P 7T P + 

xy yz 1 


7T P 7T P - 

yz xy 


7T P 7T P - 

yy xz 


7T P 7T P 

xz yy 




2 






-^xx 


7r xz 7r 0z " 


-1- 7T P 7T P 


4- 7T P 7T P 


-1- 7T P -7T P 

1 "Oy n xy 






2 




2 P = 

yx 


7T P 7T P - 

''yz^Oz 


+" 7r 0z 7I 1/z 


+ TT P y Tt 0y 


+ ^Oy^yy 






2 






^zz^Oz " 




+■ 


+" 






2 




JP _ 


1- E 




+ (^z) 2 

2 


+ 0(tt 3 ) 



0(tt 3 ), 

+ o(^ 3 ), 

+ 0(^ 3 ), 

+ o(^ 3 ), 

(3.49) 



Now, using (I3.3[) we obtain the expression of S a ,V a , lZ a b 

in terms of 7r^ 5 , which we substitute into the effective 

Hamiltonian f|2.17|) . 

In this way we derive the effective Hamiltonian of the 

Goldstone modes in terms of the canonical sets of 7r p , 
~ My 
and 7f p . In the momentum space it reads 



0,x,y ,z 



where 



w = n\ + nl + nl, (3.50) 

Hi = Cl^l yM + 3£< fc ^ fc , (3.51) 

n = A p n p l k n^ k + Bl^t >k nl !k , (3.52) 

H P = (^ p )tX p ^ p , (3.53) 
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with 



2 V 
Po 

— -k - 

Po 

2J P 
-±k 2 
Pa 



^SAS 



2^W 

AsAS 



Asas 



2^/1^2 



+ £cap(l-/3 2 ), 



jf = (i-/3 2 )j s + /3 2 j; 



p _ 



I 77" 7T 7T 7T 

V yy,ki n xz,kt "xj/,fe' yzM 

A p k A z /2 
Az/2 B£ 



The canonical commutation relations are 




(3.54) 



~P ~ p 

7T 7T 



= iS(k + k') 



(3.55) 



for each fi = 0, x, y, z. 

We first analyze the dispersions and the coherence 
lengths of the canonical sets of the modes 7Tg y and ttq z 
from (I3.51[) . Since the ground state is a squeezed co- 
herent state due to the capacitance energy e cap , it is 
more convenient [l| to use the dispersion and the coher- 
ence lengths of tTq v and tTq z separately. The dispersion 
relations are given by 



E *J* = ^k 2 
Po 

— -k z 



A 



SAS 



Asas 
Po 2 V / T 3 £ 2 ' 

and their coherence lengths are 



+ e cap (l-/3 2 ), (3.56) 



(3.57) 



Coy =2l B ^ 



-- + 2e cap (l-/3 2 )' 



'nJ*y/l -fi 2 



*SAS 



(3.58) 
(3.59) 



They describe a ppin wave. 

The similar analysis can be adopted for the canonical 
sets of n% and tt^ z in (|3.52p . The dispersion relations 
are given by 



2V 
-±k 2 

Po 

2 -^k 2 

Po 



*SAS 



A 



SAS 



-2e x (l-/3 2 ), (3.60) 



(3.61) 



Their coherence lengths are 



2l f 



TTjf 



A 



Asas 



-0 2 



4^(1 ~P 2 Y 



iSAS 



(3.62) 



(3.63) 



It appears that is ill-defined for Asas — > in (|3.60l) . 
This is not the case due to the relation (I3.65|) in the ppin 
phase, which we mention soon. 

Finally, making an analysis of the Hamiltonian (|3.53|) 
as in the case of the spin phase, we obtain the condition 
for the existence of a gapless mode, 



^SAS 



VT^W 2 



A 



SAS 



-4e x (l-/3 2 ) 



A| = 0. (3.64) 



It occurs along the ppin-canted boundary: See (5.3) and 
(5.4) in Ref . [20] . Inside the ppin phase, since we have 



^SAS 



A 



SAS 



^J 2 



4ex(l-/? 2 ) 



A| > 0, (3.65) 



there are no gapless modes. 



C. CAF phase 

Finally we analyze the CAF phase. This phase is char- 
acterized by the order parameters (|2.20p . which we may 
rewrite as 

(3.66) 



where 



C n 



COS 



sm( 



a, 



cg p = cos 9p = \J\ - /3 2 , 



ce s = cos ( 



= V 1 - " 



s 6/3 = sin 9p = -/3, 

s 9s = sinf 



A. 



SAS 



A 
(3.67) 



. The order parameter ijn) is quite complicated. Never- 
theless, the problem is just to find an appropriate rota- 
tion in the SU(4) space so that the order parameter has 
only a single nonzero component after the rotation. 

There are two ways. One is by choosing the rotational 
transformation as 



U^p = exp[i9 s T yz }exp[i9 a T xy }Vp(9p), 
with Vp given by f|3 . 38[) . and we obtain 



(3.68) 



(3.69) 
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In this rotated basis, the further analysis goes in parallel 
with that given in the spin phase. Another choice of the 
rotational transformation is given by 



U l,0 = CX P 



= exp 
obtaining 



7T 

2 y 



exp 



T pc(0) 



U p 



j.l.U 



exp 

7T 

-i—T T 
2 



r(0) 



(*--s) 



(3.70) 
(3.71) 



In this rotated basis, the further analysis goes in parallel 
with that given in the ppin phase. We call the rotated 
basis of the SU(4) group given by (|3.68|) . the s-coordinate, 
and the rotated basis given by (13.701) . the p-coordinate. 
They give the identical results. 

We make an analysis by employing the s-coordinate. 
Namely, we define the SU(4) isospin operator in the s- 
coordinate by 



exp 



75 



j«c(0) 



(3.72) 



for the fields in the ppin phase. 

We require (|3.72|) to satisfy the SU(4) algebraic rela- 
tion, 

[T x l(x,t),T y l(y,t)] = ip^ (x,t)S(x - y), (3.78) 

from which we obtain the canonical commutation rela- 
tion, 



[tt- (a,*),*- (y,*)] = iS(x-y), 



(3.79) 



with 7f^ 



We are able to derive the effective Hamiltonian for 
the Goldstone modes precisely as we did for the ppin 
phase. Namely, we obtain the relations between the orig- 
inal fields I^jy and the fields irf? s from (13.72j) . We give 
the explicit relations in Appendix: See ()B2j) . and (|B1[) . 
Thus we derive the effective Hamiltonian of the Gold- 
stone modes in terms of the canonical sets of 7r^ and ff^ . 
Working in the momentum space, the effective Hamilto- 
nian reads, 



where 



Hf + U\ 



(3.80) 



7 5 
7<5 



(3.73) where 



with ([321) and (j3~45|) . 

The eight Goldstone fields are, 



rsc(i) _ 



r sc (!) — ^ sc 



and the symmetry breaking pattern reads 
SU(4) -> SU(2) Z - ® U(l)2j. ® SU(2)x; 



(3.74) 



(3.75) 



just as in the cases of the spin/ppin phase. 

Here we remark how the Goldstone modes in the CAF 
phase are transformed into those in spin/ppin phase at 
the phase boundary. On one hand, the field ir s ^ v shift 
smoothly to the field (|3.13[) . by the inverse transforma- 
tion of (|3.68|) . or by taking the limit a, f3 — > 0, as 



(3.76) 



so that subscript of tt s ^ v perfectly matches with tt^ for 
each [iv in the spin phase. On the other hand, 
shift smoothly to (13.43|) . by the inverse transformation 
of exjp(i9$T yz )exip(i9 a T X y) , or taking the limit a — > 1 as 



1 zz ' 



zy 



7T — > 7T H 



(3.77) 



with 



— ^i,fc(^<c0,k) ^So.fe + ^2,fe(^o,fe)^yo,fe! (3.81) 
= C^r^fe, (3-82) 



and 



•) AnC fl 1 



Lxj k — — j 1 K 

Pa * 

2 o o « n - 4(Sn Cfl + Cfl )£y 

( - z 2,fc — l c e a J s t s e s J i ) K "I ^ i 

Po ^ 

J? = c 2 9a J s + s 2 8a Jf, M = Ac 2 9a e x + Aoc, 1 , (3.83) 



xz,k 

VV,k 
-sc 
yx.k 



M% c = 



( A c 
c c 
e c 



V 



& 

C c 

r 






F c 



B c 








a 
If 







a c 
D c 
d c 



\ 




6 C 
d c 
E c ) 

(3.84) 
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The Matrix elements in (I3.84[) are given by 



D. CAF phase in Asas 



A c = 

B c = 

C c = 

D c = 

E c = 
+ 

F c = 
and 



2 T /3 , 2 jd 



,2 T/3 j „2 t/3 



c e a J 3 



*9c, J l 



M 

T ~ 

Ap 



Jf + ^ -2cjLe 



2fc2 

Po 

Po 

2*^ 

Po 
Po 
Po 

222 i)/ 2 2 i 2 \ 2 

s 0p S s C 8 a e c&p ~ z { c Bp S 8 s +Cf,,)S a e 



zs e /3 c e s e X' 



M 
T 



The effective Hamiltonian in the CAF phase is too 
complicated to make a further analysis. We take the 
limit Asas — > to examine if some simplified formulas 
are obtained. In particular we would like to seek for gap- 
less modes. Such gapless modes will play an important 
role to drive the interlayer coherence in the CAF phase. 

In this limit we have 



2 / 2 t/3 , 2 jP\ , 2 jc 

C 9 S \ s e a J 3 < c e a J i J ' s e s J i 

2 (2 jp , 2 jP\ , 2 tc 



Ap 

2c Sfi 
M 

Y 



C 8 s S 8{, e a 



c :1 



-*SAS 

A z : 



±W1 



-ASAS 

A z 



Po 



X' 



71/ 



(3.85) 



c<? 5 = c 6a , s 8s = se Q , Aoc^ 1 = A z 
a c = 6 C = c c = e c = L = 0. 

By using the above equations, (13.81)) become 

4 



— Jffe 2 + A z 
Po 



2fc 2 



Po 
2fc 2 



-C8 g C29 a J2 



S 28 B C 8 S 



with 



S0 5 S20„J 2 -r-^-r 

Po 

2fc 2 T0 ^ 

ce s J2 +s 2 e ce s e x 

Po 

S28 a S26 s ,2k 2 



iSAS 



4A f 



-Ce a S26g^a, 



Vi.k 



tO,k 



n sct T7 sc 
J /l,fe ? 'l,fe' 



'yo.fc 



V2 



(3.88) 



(3.89) 



(3.90) 



Po V 



t/ e (/q Js 



From p.89p we have the dispersion and the coherence 
length for mode 77^° 



2 N 

S 8f,\^ e X ~ £ cap)J „""' 



£"1 



Jffe 2 + A z , £ 



iV 



Po 



'/I 



2, L! J^l. (3.91) 



(3.86) 



with 



This mode is reminiscent of the spin wave (|3.24p in the 
spin phase. 

We next investigate Uf in f|3.82|) . It yields 



t/3 _ 2 t , 2 T d t/3 _ S2fl g f jd J \ 
J l — c 8,j J s T Sg^^s i J 2 — 2 ' s s '' 

t/3 _ 2 xd , 2 j ja _ 2 jd , 2 t 



L = — 



s 26 b 



£ cap) + C0 ( 



A 



SAS 



A 





_ njsc 
— H 2,l 


-T/SC 

rt 2,l 




' 4 

.Po' 


n 2,2 


-»sct 

= *2,i- 



-Jffc 2 + A z 



scf sc 



N = 



S28 s S28 a S g(j _ 

4ci_ey + 2sfl_e c 



Asas , _ 2 



{c0 5 ce a s 6l3 e a + A z ) , where 
(3.87) 



vT.k 



yx.k 



where we denote S29 a = sin2# Q , 529,3 = sin 26 p, and 
528,5 = sin 29$. 

It can be verified that the effective Hamiltonian (13.81)) 
and p.82j) reproduce the effective Hamiltonian in the spin 
phase p. 201) . by taking the limit a — > first, and then di- 
agonalize this Hamiltonian with the transformation Vg 1 , 
or taking a, j3 — > 0. On the other hand, we reproduce the 
effective Hamiltonian in the ppin phase (|3.50p , by taking 
the limit a 1, in (|3~5Tj) and (|3~Mj) . 



V2 



(3.92) 
(3.93) 
(3.94) 

(3.95) 



and 



l 2,fc 



/ K%k \ 

Z.SC 

yy,k 

~sc 

yz,k 

V ^xy,k / 



M 



2.2 



/ c c r \ 

f c F c 

D c d c 

\ d c D c J 



(3.96) 
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with 



where 



2fc 2 
Po 
2k 2 



C c = —Jt + ^ + 2e x (cl ^as 



Po 
2k 2 



Az 
2 

A* 
2 



A 2 



d c = —4e a Us - Jf) - ^ + ^«e ca p - 2(1 + s 2 Je x ] 
Po ^ 



f c _ Az 2 / A|a; 



(3.97) 



4fr 2 

-E^ = jf + 2A z + 8c 2 e x , 

Po 



E vl° = | fc | / M (^(c^ J* + s 2 e J s ) + 2s 2 en (e D e x ) 
V Po V Po 



(3.104) 



The annihilation operators ?^ c fc (i = 3, 4) are given by 



^xy.k + ^rrz.fc' 



From p.93[) we have the dispersion and the coherence 
length for the mode jyjf, 



= —J?k 2 + A z , f>5 c = 2l B J^, (3.98) 
Po V A z 



which have exactly the same value as (|3.91[) . 

We next analyze %| C 2 and take Asas = for the sake 
of the simplicity. This Hamiltonian can be diagonalized 

as 



/ Klk \ f ( AS \ 



T/SC 

rt 2,2 — 



fr sc 
K yz,k 
~ sc 



X S y y 

A* c 2 
V OA-/ 



7T SC 

yy.k 

K yz,k 
- sc 
\ W xy,k / 



(3.99) 



where 



2tf_ 

Po 
2k 2 



J s d + A z +4c 2 Q e x , 



\sc . rpc jd 

yy ~ J ~ - s 



X S y z =D C +d C 



po 
2fc 2 ( 

Po 
2fc 2 



A!i = # c - <f = J„ d + A z + 4c 2 , e 



and 



Po 



a c X' 



(3.100) 



-sc _ ^zz.fc "j~ ^j/y.fc -sc 
"zz.fc ~ ^= ' "< 



yz,fc 



V2 



^yz.fc + ^xyM ^ sc 



V2 



VV,k 
xy.k 



'^xz.k + ^yy^fc 



V2 



~ n yz,k + ''"xy.fc 



V2 



(3.101) 

The fields (|3.101[) satisfy the commutation relation 

K v ,k, = + k'), [n y c z k , 7^,] = iS(k + k'). 

(3.102) 

We can rewrite the Hamiltonian (|3.99p as 



Hf 2 = / d 2 kV. s 2 % = I 



set sc 



(3.103) 



", c „sct „sc 



Vi.k 



V2 \ ]j \y V \z 

They satisfy the commutation relation, 



T7 SC T7 SCt 



= 5 ij 5{k-k'), 



(3.105) 



(3.106) 



with i, j = 3, 4. 

We summarize the Goldstone modes in the CAF phase 
in the limit Asas —> 0. It is to be emphasized that there 
emerges one gapless mode, ?74° fc , reflecting the realization 
of an exact and its spontaneous breaking of a U(l) part 
of the SU(4) rotational symmetry. Furthermore, it has 
the linear dispersion relation as in (13.1041) . as leads to 
a superfluidity associated with this gapless mode. All 
other modes have gaps. 

We comment on the existence of the two modes (|3.90[) 
and (I3.95[) . Their dispersions (|3.91[) and (|3.98[) are similar 
to that of the spin wave (|3.25[) . The difference between 
the dispersion of these two modes and the spin wave is 
the stiffness dependence. (|3.9ip and p.98[) have the stiff- 
ness structure of the linear combination of the intralayer 
stiffness and intcrlaycr stiffness. This can be understand 
because the CAF phase has a layer correlation. The other 
modes are massive due to the Coulomb energy and the 
Zeeman gap. 



IV. DISCUSSION 

We have presented a systematic method based on the 
formula (|3.2[) to investigate the symmetry breaking pat- 
tern and to derive the effective Hamiltonian for the Gold- 
stone modes in the v=2 bilayer QH system. There are 
eight Goldstone modes in each phase, which are shown 
to be smoothly transformed one to another across the 
phase boundary. In particular, we have analyzed the 
CAF phase in detail. 

The interlayer phase coherence and the Josephson ef- 
fect are among the most intriguing phenomena in the 
v = 1 bilayer QH system [lj. They are enhanced in the 
limit AgAs 0. It is natural to seek for similar phenom- 
ena in the v — 2 bilayer QH system. We may naively 
expect them to occur in the ppin phase. However, as we 
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have found, almost all electrons are moved to one of the 
layers in this limit. 

This is not the case in the CAF phase, where the elec- 
tron densities can be controlled arbitrarily in both layers. 
In the CAF phase we have investigated the dispersion re- 
lations and the coherence length in the limit Asas — > 0. 
Remarkably, we have found one coherent mode whose 
coherence length diverges. Furthermore it has the lin- 
ear dispersion relation. It might be responsible to the 
interlayer phase coherence. 



where I2 is the unit matrix in two dimensions. Nine 
remaining matrices are simple products of the spin and 
pseudospin matrices: 



T spin r ppin 



7 -spm T ppin 



T a 

T a 

T a 

-T a 



spin ppin 
a 'y 



-~lT a 
%T a 



(A5) 
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We denote them T a0 = ± T f la , T 0a = §7-^, T ab = 
h T a Pin ' T f Pm ■ They satisfy the normalization condition 



and the commutation relations 



(A6) 



(A7) 



Appendix A 



The special unitary group SU(N) has (N 2 — 1) gen- 
erators. According to the standard notation from el- 
ementary particle physics [28], we denote them as Xa, 
A = 1,2,..., N 2 — 1, which arc represented by Hermi- 
tian, traceless, N x N matrices, and normalize them as 



Tr(AAA B ) = 2S A b- 
They are characterized by 

[Xa, As] = 2ifABcXc, 
4 

{Aa,A_b} = —2dABcXc, 



(Al) 



(A2) 



where Jabc and (Iabc are the structure constant of 
SU(N). We have A^ = ta (the Pauli matrix) with 
Iabc = ^abc and (Iabc = in the case of SU(2). 

This standard representation is not convenient for our 
purpose because the spin group is SU(2) x SU(2) in the 
bilayer electron system with the four-component electron 
field as * = (V> f ^ fJ -,V> bt ,^ H )- Embedding SU(2) x 
SU(2) into SU(4) we define the spin matrix by 



spin _ 



Ta 
T a 



(A3) 



where a = x,y,z, and the pseudospin matrices by, 



ppin _ 



l a 

1 2 

1 2 

-1 2 



ppin _ 







-iU 




where fuvjS aiu 1 is the SU(4) structure constants in the 
basis (IA3|) - (IA5|) . Greek indices run over 0, x, y, z. 



Appendix B 



We express the rotated isospin fields 23, in terms of 
the eight Goldstone fields t s ^ v up to the second order, 



n~SC 

X 02 — 



'7"SC _ 



JSC 



sc sc , 
xz"yx 1 


n sc sc _ 
1 yx xz 


7T SC 7T SC - 

yz ' xx 


7T SC 7T SC 

xx yz 

— 4 




-0(7T 3 ), 
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sc sc , sc sc _ 
xx yy 1 yy xx 


n sc sc _ 
xy yx 


7T SC 7T SC 
yx xy 


-0(t 3 ), 




2 






_sc _sc 

y x yO 


4- ir sc -7r sc 
' ~yO~yx 


1 ^-sc sc 
1 n xx n xO 


T n xO n xx 


+ 0(n 3 ), 






2 




7T SC 7T SC 

yy yO 


1 sc sc 

"f KyQTlyy 


1 sc sc 
T ^xy^xO ' 


_i_ ,_sc sc 
' n x0 n xy 


+ 0(n 3 ), 






2 




_sc _sc 
_ ^yz^yO 




1 _sc sc 


4- TT SC 7T SC 


+ 0(tt 3 ), 






2 




sc sc , sc sc _ 
"xy yz 1 yz n xy 


^sc sc _ 
yy xz 


7T SC 7T SC 
n xz"yy | 


■0(n% 




2 







2S = 1 _ j2 [ ^ )+ [ - ^ +0 {**). (Bi) 



fi=0,x,y,z 



(A4) We next give the relation between the original isospin 
field I^v and the rotated field I s £ v in the s-coordinate of 
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the CAF phase. with (f3T67|) . 

Iqx = c 8a c gfl ce s T c x + c^s^Z^ - se a sg f) l s xx + sg a c gfj cg 5 X x c z 

1-Oy = Ce s ^Qy + Sg s Xy C x , 

Zoz = - c e a s e/1 Cg 5 2Q x + Cg a Cg e T?Q z - Sg a c e/1 l xx - Se a Se fi Cg s T x c z 
ZrO = Cg s T s x0 + Sg 6 Tf z , 

Zxx = -Se«Sg p Cg 6 Tl x + Sg a Cg m + Cg a Cg p l s xx + Cg a S g ^Cg s l x ° z 

+ s B a Sg f) Sg 5 I y c y + Cg a Sg^Sg s l s z0 , 

J — T sc 
-L-xy — -L-xyi 

Ixz = -~Sg a C ei3 Cg s lQ x - S^S^Zqz - Cg a S ei3 l xx + Cg^Cg^Cg^l^ 

i TSC I 'T-SC 

+ s o a C0 S0 6 x yy + cg a ce !3 se s x z0 , 

yO — Cg a -L yQ - Sg a ± zyl 
Zyx = -COpSg^ly + Cg p Cg 5 l s yx + Sg p l s yz , 
X yy = c e a Se s I^ C x + Sg a Sg 5 l s xz + Cg a Cg s l y c y - Sg a Cg 6 l s zQ , 
lyz = S 9f) Sg 5 m - S 6l3 Cg s l y C x + Cg^lyl, 

zO = Sg a Sg s l 0x - Cg a Sg 6 l xz + Sg a Cg s -L yy + Cg a Cg 5 ± zQ , 
Izx = -S S? Sg s I^Q + Cgpl S zx + Sg^Cg^H, 

T-zy = s e a Z y o + c e a Zly, 

T-zz = -cgpSg s 2 x c - s efj l s z c x + c e pCg 5 l s z c z , (B2) 
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